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Abstract 

Neutrinoless double beta ((3/3q u ) decay is considered within bilinear R-parity 
breaking super symmetry, including the full one-loop corrections to the neutrino- 
neutralino mass matrix. Expected rates for (3f3o u decay in this model are discussed 
in light of recent atmospheric and solar neutrino data. We conclude that (a) tree- 
level calculations for (3(3q v decay within the bilinear model are not reliable in the 
range of parameters preferred by current solar and atmospheric neutrino problems. 
And (b) if the solar and atmospheric neutrino problems are to be solved within 
bilinear R-parity violating SUSY the expected rates for @Po u decay are very low; 
the effective Majorana neutrino mass at most 0.01 eV and typical values being one 
order of magnitude lower. Observing P(3q u decay in the next round of experiments 
therefore would rule out the bilinear R-parity violating supersymmetric model as 
an explanation for solar and atmospheric neutrino oscillations, as well as any hier- 
archical scheme for neutrino masses, unless new neutrino interactions are present. 
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1 Introduction 



Neutrino physics has entered a new era recently with the announcement by the Super- 
Kamiokande collaboration of rather conclusive evidence for neutrino oscillations [1] in 
atmospheric neutrino measurements. This experiment, together with the oscillation in- 
terpretation of the long-standing solar neutrino puzzle [2] now provides important infor- 
mation on neutrino masses and mixings and may-be the first look to physics beyond the 
standard model [3,4]. 

However, neutrino oscillation experiments, while being extremely valuable, can not answer 
two fundamental questions in neutrino physics. First, they are only sensitive to mass 
squared differences and thus can not fix the overall mass scale of neutrinos. And, second, 
due to the V-A nature of the weak interaction neutrino oscillations can not distinguish in 
practice between Dirac and Majorana neutrinos. Other experiments on neutrino masses 
are needed in order to reconstruct the neutrino mass matrix. Neutrinoless double beta 
decay is a prominent example of such kind of experiments. 

Neutrinoless double beta (PPo v ) decay has for a long time been known as a sensitive 
probe for physics beyond the standard model (SM). Non-observation of (3(3 Qv decay has 
been used to derive stringent limits on various extensions of the SM, like, for example, 
left-right symmetric models [6], leptoquarks [7] and supersymmetry [8,9,11]. However, 
PPov decay has yet to be observed experimentally. 

Although there might exist a variety of mechanisms inducing PPq v decay in gauge theories, 
one can show that whatever the leading mechanism is at least one of the neutrinos will be 
a Majorana particle [12]. The observable in PPq u decay, the effective Majorana neutrino 
mass, is in general a superposition of different mass eigenstates: 



where U e j characterizes the couplings of the mass-eigenstate neutrinos to the electron in 
the charged current and the prime indicates that the sum runs over light mass eigenstates 
only. If neutrinos have non-zero mass, also non-zero mixing among them has to be ex- 
pected, so that in general (m v ) does not coincide with the electron neutrino mass probed 
in tritium beta decay. 

Currently the most stringent experimental bound [13] gives an upper limit of the order of 
(m„) < 0(0.2 — 0.5) eV. There exist two independent proposals for future experiments 
which might improve the sensitivity on {m v ) by up to one order of magnitude or more 



Here, we concentrate on the calculation of expected rates for PPq u decay within bilinear 
R-parity violating (BRPV) SUSY. While PP 0u decay has already been considered in the 

4 The oscillations which are Dirac-Maj or ana-sensitive must violate lepton number by two units 
and are helicity suppressed [5] 



( m u) = U ej m 




J 



[14,15]. 
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literature before within the explicit BRPV SUSY model [9-11], it has so-far only been 
treated in lowest order of perturbation theory considering the neutrino-neutralino mass 
matrix only at the tree-level approximation. Here, we take into account the full one-loop 
corrections to the neutrino-neutralino mass matrix and especially concentrate on those 
regions in parameter space in which the model can solve simultaneously the solar and 
atmospheric neutrino problems [16]. 

We have found that there exist important regions in the parameter space of the model - 
namely those where the BRPV SUSY model can account for the solar neutrino anomaly 
through matter-enhanced oscillations - where the tree-level estimates for f3/3 0v decay 
fail rather badly. Thus the one-loop corrections considered here play an important role 
in BRPV SUSY. Their inclusion is definitely necessary in order to predict reliably the 
effective Majorana neutrino mass relevant for /3/3 Qu decay in a way consistent with the 
results from present oscillation experiments. 

This paper is organized as follows. In the next section we set up the notations and discuss 
the model at tree-level. Then, we outline briefly the extension of the calculation including 
the one-loop corrections. Further details for these can be found in [17]. Section 4 discusses 
our numerical results. 



2 Bilinear R-parity violation and neutrino mass at tree-level 



In the following we use conventions such that in the limit were the R-parity violating 
parameters vanish the usual MSSM notations of refs. [18] are recovered. For the BRPV 
case see ref. [19,20] for the conventions we adopt. The supersymmetric Lagrangian is 
specified by the superpotential W given by 



W = e ab WiQ^M + h%Q\D^l + h%L\R 3 H a d - ^H a d H b u + e^H* 



(2) 



where i,j = 1,2,3 are generation indices, a,b = 1,2 are SU(2) indices, and e is a com- 
pletely antisymmetric 2x2 matrix, with E\2 = 1. The symbol "hat" over each letter 
indicates a superfield, with Qi, Li, H^, and H u being SU{2) doublets with hypercharges 
|, —1, —1, and 1 respectively, and U, D, and R being SU{2) singlets with hypercharges 
— 3, |, and 2 respectively. The couplings hu, hp and Ue are 3x3 Yukawa matrices, and 
/i and ti are parameters with units of mass. The last term in eq. (2) is the only .R-parity 
violating term. 

Supersymmetry breaking is parameterized with a set of soft supersymmetry breaking 
terms, 



V soft = MfQTQ) + MfUV] + MfDiD* 
+m 2 H HTH a d +m 2 Hu H a u *K- 



\M S \ S \ S + \M\\ + iM'A'V + h. 



lid " " 11 u u IZ " " " Z Z 

A%Q1XJ 3 H h u + A%Q h fi 3 H a d + A%L\RjH d - BfiH a d H b u + B^H* 



(3) 
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and again, the last term in eq. (3) is the only R-parity violating term. The bilinear term 
in (3) leads in the neutral part of the scalar potential to terms linear in the sneutrino 
fields. Thus, in general the sneutrino fields acquire VeVs. This in turn leads to mixing 
between the gaugino and lepton as well as to mixing between the scalar leptons and the 
Higgs fields [20,21]. 

For our purposes the most important aspect is the neutrino-neutralino mixing, since it 
leads at tree-level to one massive neutrino state. In the basis, ^/' T = (ip^ , ip\ 2 , ty\ , 
—iXsj^jjj^ff) the neutrino-neutralino mass matrix at tree-level can be written as: 



Mo 




(4) 



Here, the sub-matrix m contains entries from the bilinear ljl p parameters, 



m 



^ -\g'v e \gv e e e ^ 
K ~\g'v T \gv T e Ty 



(5) 



where Vi :- 



(i>i) and -M x o is the MSSM neutralino mass matrix, given by, 



M x o 



Mi 





-\g'v d 


\g'v u 





M 2 


\gv d 


-\gvr 


-\g'vd 


\gvd 





-n 



\ \g'vu -\gv u -fi 



(6) 



There are two interesting aspects concerning M. . First, M.® has such a texture that at 
tree-level only one neutrino gets a non-zero mass [22], leaving two massless (but mixed) 
states in the spectrum. And second, at tree-level the neutrino mass is strictly proportional 
to the "alignment vector" |A| 2 , where, 

A := ev d + vfx. (7) 



Thus, at tree-level the individual e« and Vi are not constrained neither by the neutrino 
mass measurements nor by neutrinoless double beta decay, as long as they are sufficiently 
aligned. However, we would like to stress (more details below) that this is a pure tree-level 
result. Once the calculation is improved to one-loop order current experimental hints on 
solar and atmospheric neutrino oscillations provide rather stringent constraints not only 

— * 

on A, but also on the individual BRPV parameters, q and v j. 
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Assuming that m <C M.q one can find [9,11] a simple formula relating the effective Majo- 
rana neutrino mass to the supersymmetric parameters: 

It has been shown in [11] that within BRPV the contribution from (m u ) as given above 
is the dominant source for /3/3o v decay. In the following we will concentrate on this BRPV 
mass mechanism only, improving it by taking into account the one-loop corrections to the 
neutrino-neutralino mass matrix. 



3 One-loop corrections to the neutrino-neutralino mass matrix 

As we have seen the effective neutrino mass matrix has a projective structure, such that 
only one neutrino gets a mass at tree-level. As a result for a realistic description of the 
neutrino spectrum one has to improve the calculation to 1-loop order. Q A shortened 
description is given below, for a complete listing of all necessary couplings etc. see ref. [17]. 
However, most important for the understanding of the importance of the loops is the fact 
that these contributions explicity break the projectivity of the tree-level mass matrix, 
incorporating contributions which are proportional to the e, themselves, as we will show 
explicitly below. In contrast, as discussed above, the tree-level mass matrix is sensitive 
only to A. 

The full neutrino-neutralino mass matrix including the 1-loop corrections is given by 

Ma = Mf e + AMy, (9) 

where AMy are the 1-loop corrections defined by 

AMy = i(n«(p?) + n^ 2 ) - m^M) - m^aipj)) (10) 

where E y - and n y - are self-energies. There are three simple topologies of relevant Feynman 
diagrams contributing to the neutrino-neutralino mass matrix [17]. 0. Here, DR signifies 
the minimal dimensional reduction subtraction scheme and is the renormalization 
scale. As pointed out in [17] the inclusion of the tadpole diagram is essential in order to 
obtain gauge invariance of the calculation. 

Figure 1 shows the relevant Feynman graphs. Internal particles in the scalar self-energies 
can be either (q — q), (charged scalars-charginos) or (neutral scalars-neutralinos) , for the 

5 With two massless neutrinos, one angle of the neutrino sector of the theory could be rotated 
away. Thus a discussion of the predictions of the theory for the solar angle is meaningless at 
tree-level. 

6 For a complete description see ref. [17] 
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Fig. 1. Simple one-loop topologies contributing to the neutralino-neutrino mass matrix, see text. 

gauge loops it can either be (W ± — charginos) or (Z°— neutralinos). Which of the loops is 
most important depends both on parameters and whether one considers the heavy states 
("neutralinos") or the light states ("neutrinos"). Here we concentrate on the "neutrino" 
states. For these only the (d — d), (charged scalars-charginos) and (W ± — charginos) com- 
binations do indeed contribute. For large values of tan (3 generally the (d — d) loops are 
most important. We will therefore concentrate on this loop in the following, noting in 
passing that the basic structure of all the self-energies are the same and can be found by 
replacing internal masses and couplings correspondingly [17]. 

It is interesting to note that the tree-level result of neutrino masses being strictly propor- 
tional to |A| 2 is no longer valid once the one-loop contributions are taken into account. 
This can be shown for example for the down-type squark loops, for which H^pf) and 
Eij(p-) are given by, 

KiM) = ^ E {Ol%O d L,L + O n R % s Otkis)rn k B (ml ml m 2 s ) (11) 

k,s 

^M) = ^ E {° n R%sOtL + O n L % s O d ^is)Bi{ml ml ml) (12) 

k,s 



where B and B 1 are Passarino-Veltman functions [24], m k and m s are the down- type 
quark, down-type squark masses and the various O are neutralino-quark-squark couplings, 
in our notation given by, 

R T jk = (^)We " I tan^A/^Rf^Rl^ - (h^R^R^^ (14) 

where the h d denote the down-type Yukawa couplings and 0^ k = (0 R ^ ik j and R d ^ k = 

i^Ljik) ■ The rotation matrices R d and R^ are the ones which diagonalize the quark and 
squark mass matrices, respectively, while Af diagonalizes the neutralinos/neutrinos. 

That terms proportional to e$ survive in eq. (11) is most easily seen assuming the BRPV 
parameters are small, as suggested by the present indications from solar and atmospheric 
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neutrino data. Then one can block-diagonalize the neutrino-neutralino mass matrix per- 
turbatively at tree level in terms of the expansion parameter £ = m ■ -M~o [23] as, 



V T (l - Ifitf) -V T £ 
Af*=[ ^ | (15) 

where N* is the matrix diagonalizing the MSSM part of the neutralino mass matrix and 
Vj describes the mixing of neutrinos among themselves. 

The full form for the expansion matrix £ can be found, for example, in [11]. For our 
purposes it suffices to state that in the limit A = the matrix Vj is diagonal, and all 
elements of £ vanish except £ i3 , which take the simple form, 

& = -- (16) 



Inserting this result for A = and for simplicity considering only i,j = 1,2,3, (Ejj 
vanishes for i,j = 1, 2, 3 in this limit) 11^ can be written as, 

KM) = ^^E( R f, fc+ 3Ri + ^0l(M^l^oK 2 ,m 2 fc ,m2), (17) 

i07r ^ k,s 



where, for simplicity, we have assumed that hd is diagonal. Eq. (17) demonstrates that 
the entries in in the "neutrino sector" are proportional to e^e^. This shows explicitly 
that in the limit where the tree-level neutrino mass vanishes the loop contributions do 
not and can, potentially, be rather important. Moreover, from this example we can draw 
two conclusions. First, 1-loop contributions break the projectivity of the mass matrix 
{mfj ee ~ AjAj at tree-level) and thus the degeneracy of the two lightest states is lifted. 
And, second, the size of the ratio of the 1-loop to the tree-level entries of the mass matrix 
should be controled mainly by the quantity |e| 2 /|A|. 7 



4 Numerical results 



In our numerical study we assume unification at a scale Q = My with standard minimal 
supergravity boundary conditions, 



A f = A h = A, = A 



7 In the numerical calculation we have found that this is indeed the case. However, the loops 
depend also strongly on tan f3, because large tan (5 leads to large Yukawa couplings in the down 
sector, and as shown in eq. (17) the 1-loop entries strongly depend on hd- Numerically, variations 
of other SUSY parameters have been found to be much less important. 
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Fig. 2. Effective Majorana neutrino mass as a function of Amj 2 for data points which have 
s'm 2 (28 so i) > 0.6 and solve the atmospheric neutrino problem. 



M 3 M 2 M : Mi/2 ■ 

We run the RGE's from the unification scale Mjj ~ 2 x 10 16 GeV down to the weak 
scale, giving random values to the fundamental parameters at the unification scale. We 
then check that the numerical values obtained from the RGE running correctly break 
electroweak symmetry. Moreover, we accept only those points for further study, which 
fulfill phenomenological constraints from negative Higgs and SUSY particle searches at 
accelerators [25]. 

Although this procedure is not essential for the calculation of the neutrino masses in the 
model, it allows us to reduce the number of free parameters considerably and can be 
viewed as a test for self-consistency of the parameter ranges under consideration. 

For the Jji v parameters, we use the constraints from solar and atmospheric neutrinos found 
in [16,17]. These two sets of measurements imply that BRPV parameters have to be small, 
i.e. |e| and |A| should be smaller than O(GeV) and O(0.2GeV 2 ) respectively for typical 
MSSM parameters smaller than, say 1 TeV. Moreover, measurements of (or limits on) 
neutrino angles fix (or yield limits) on ratios of R-parity breaking parameters. Here 
we summarize these restrictions as follows [17]. The atmospheric neutrino measurements 
require A M ~ A r , whereas the negative results from the CHOOZ [26] and Palo Verde 
reactor [27] experiments require that A e should be smaller than A e < (9(0.3) ^A^ + A^. 
The solar neutrino problem can be either solved with relatively large mixing (LMA-MSW 
or vacuum oscillations), which implies that all q should be similar, or by small mixing 

8 Although smaller than usual supersymmetric parameters, such a suppression might be actually 
expected in scenarios with radiative R-parity breaking ref. [17,19] 



B = Bi = A - 1 , 





M 2 Qi = M* t = M% = ml , 
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Fig. 3. Ratio of 1-loop corrected effective Major ana neutrino mass to its tree- level value as a 
function of Am 2 2 for data points which have sin 2 (2# so ;) > 0.6 and solve the atmospheric neutrino 
problem. 

(the SMA-MSW solution), the latter implying e e ~ (few)l(T% )T . 

We have determined the expected values of {m v ) as a function of Am 2 2 for about 10 4 
calculated points, which solve the atmospheric neutrino problem. Predicted values of 
(m„) are rather small, reaching at most 10 -2 [eV] for the large mixing solution (LA- 
MSW) of the solar neutrino problem, as can be seen from Fig. (2). For the case of vacuum 
oscillations (ra„) will be even much smaller, around 10~ 4 [eV], as seen from the figure. 

Let us now discuss the crucial importance of the loop corrections to the neutrino masses 
in this context. In order to do this we have calculated ratios of (m v ) including the 1-loop 
corrections divided by its tree-level value. In figure Fig. (3) we show our results. As can be 
seen, if Am 2 2 lies in the range required for vacuum (or just-so) oscillations the tree-level 
and the 1-loop improved (m v ) are rather similar, whereas for larger Am 2 2 in the MSW 
range one has a substantial change from the tree-level result. Thus, tree-level calculations 
of (m v ) are certainly not accurate in this case, and the 1-loop corrections considered here 
play an essential role. 

Let us now analyze the remaining oscillation possibility to solve the solar neutrino prob- 
lem, namely the small-angle MSW solution. In this case one finds a suppression in the 
f3f3o v rate, as can be seen in Fig. (4). This result is easy to understand conceptually, as 
the /3/3q v rate must be given in terms of the only L e violating parameters in the model A e 
and e e , while sin 2 (29 so i) — > as A e , e e — ► 0. 

To close this section we mention that, although we have worked within the framework of 
a concrete model in which ft p constitutes the origin for neutrino mass and mixing, our 
conclusions are more general. In fact the smallness of effective Majorana neutrino mass 
(m„) holds in any hierarchical model of neutrino mass, of which our bilinear ]ji p breaking 
model is a particular case. Note, that although it is possible in the BRpV model to have 
two neutrinos nearly degenerate once the 1-loop contributions are included, it is never 
possible to have all three neutrinos degenerate [17]. Moreover, such points are extremely 
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1CT 3 1CT 2 1CT 1 

sin 2 (20 soi ) 

Fig. 4. Expected 1-loop corrected effective Majorana neutrino mass (m v ) as a function of 
sin 2 (2# so ;) for those points which solve the atmospheric neutrino problem. 

rare in parameter space and not protected by any symmetry in our model. In hierarchical 
models, however one expects that the maximum allowed value of (m u ) (which is achieved 
for the LA-MSW solution) can be estimated by: 

(m v ) = JX-"*i ( 19 ) 

i 

~UyAml l + UyAml tm < -J^eV 2 + 0.05yJlO-*eV 2 ~ O.OleV, 

which our numerical results confirm for the BRpV model explicitly. Note, that eq. (19) 
gives us only an upper bound on (m u ), but no lower bound and no prediction for (m u ). 

One interesting way to avoid this upper bound is the possibility of neutrinos being closely 
degenerate in mass. According to our results, this would be a clear indication that BRpV is 
not the underlying mechanism for generating the solar and atmospheric neutrino masses. 
Another is if other more exotic mechanisms for solving the neutrino anomalies are enter- 
tained, such as flavour changing interactions or decays [28]. 



5 Summary 

We have calculated the one-loop corrections to the /3/3 0l , decay observable (m u ) in bilinear 
R-parity violating supersymmetry, following the procedure developed in [17]. Since it has 
been shown in [16,17] that the model is able to solve the solar and atmospheric neutrino 
problems under certain, relatively simple assumptions, special emphasis has been put in 
our analysis on those "successful" regions of parameter space. 

There are two main results of this study. First, one- loop corrections are important for 
estimating /3/3q v decay rates in bilinear BRPV SUSY. This is due to the fact that the 
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model at tree-level has two massless states in the spectrum. This degeneracy is lifted 
once the one- loop corrections are taken into account. Since tree-level and one- loop masses 
depend on different combinations of BRPV parameters, which are a priori unknown, the 
loop corrections can be easily as big as the tree level masses. Especially this is true in those 
parameter ranges, where the model is able to solve the solar and atmospheric neutrino 
problems. 

Moreover we show that, if bilinear R-parity violating is indeed the solution to the solar 
and atmospheric neutrino problems, than the expected values of (m u ) are very small, 
certainly smaller than 10~ 2 eV, and probably even smaller than 10~ 3 eV. 

Although this conclusion might appear rather discouraging for the experimentalists, we 
would like to stress that, on the other hand, discovering (3p 0u decay at a level significantly 
larger than (m„) = 10~ 2 eV would be sufficient to rule out our model as an explanation 
for the atmospheric and solar neutrino problems. This conclusion also applies to any 
hierarchical scheme for neutrino masses. The only possible way this conclusion might be 
evaded is to consider the presence of exotic neutrino properties, such as flavour changing 
interactions or decays [28]. 
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